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Abstract 

It is very important to know the strength of annihilation contribution in B charmless 
nonleptonic decays. Bg vr^vr^ process could serve a good probe of the strength. We 
have studied the process in QCD factorization framework. Using a gluon mass scale 
indicted by the studies of infrared behavior of gluon propagators to avoid enhancements 
in the soft end point regions, we find that the CP averaged branching ratio is about 
1.24 X 10~^, the direct CP asymmetry Ctttt is about -0.05, while the mixing-induced CP 
asymmetry quite large with the value 57r7r=0.18. The process could be measured at LHC- 
b experiments in the near future and would deepen our understanding of dynamics of B 
charmless decays. 
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1 Introduction 



Recent years many efforts have been made to understand charmless decays of B mesons, which 
provide good grounds to get deep insights into the flavor structure of the Standard Model (SM), 
the origin of CP violation, the dynamics of hadronic decays, and to search for any signals of 
new physics beyond the SM. Up to now, BaBar(SLAC) P and Belle(KEK) have already 
accumulated large set of data and have made plenty of exciting measurements. Moreover, in the 
near future LHC-b experiment, the expected number of bb events produced per year is about 
10^^, it is noted that 10% of the events would fragment to Bs mesons. This high statistics will 
allow studies of rare Bg decay modes, which will provide very sensitive tests of theories for B 
decays, electro-weak interaction models and so on. 

For nonleptonic B meson decays, the most difficult aspect lies in the computation of matrix 
elements of the effective four-quark operators between hadron states. To deal with this, a 
simple and widely used approach is the so-called factorization approach(FA) [3]. In the past 
few years, new approaches, such as the QCD factorization(QCDF)|3] and perturbation QCD 
(pQCD) schemejS] have been proposed to improve the FA on QCD grounds. 

In the most cases of B meson nonleptonic decays, the annihilation contribution carries weak 
and strong phases different from that provided by the tree or penguin amplitudes, which is very 
important for studying CP-violating observables. Meanwhile, the calculation of annihilation 
contributions is interesting by itself, since it can help us to understand the low energy QCD 
dynamics and the viability of the theoretical approaches. As argued in the annihilation 
amplitude is formally power suppressed by order AqcD/f^b in QCDF. However, the annihilation 
contribution may not be small. In a recent systematic calculation of B decaysjH], it is shown 
that the annihilation contributions could cause considerable uncertainties in their theoretical 
predictions, where the contributions are parameterized in term of the divergent integral Jq ^ ~^ 
Xa = (1 + ^?Ae*'^)ln^. In this paper, we argue that the strength of annihilation could be 
probed by measuring the interesting decay mode Bg tt+tt", which is a pure annihilation 
process. In our calculation of the scattering kernel, we will use Cornwall [7] prescription of 
gluon propagator with a dynamical mass to avoid enhancements in the soft end point. It 
is very interesting to note that recent theoretical |^ and phenomenological [H] studies are now 
accumulating supports for softer infrared behavior for gluon propagator. Besides serving a 
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probe for the annihilation, the decay has some interesting features: sizable CP violation due 
to both tree and penguin operators contributing, clear experimental signatures due to its two 
charge final states. Moreover, if few percentage of final pions are mis-identified to be muons, it 
would bring considerable uncertainties to the measurement of Bg — > fi~^fi~ at LHCb. Therefor, 
the decay deserves our theoretical studies using different approaches. 

We have found that the CP averaged branching ratio of Bg tt+tt^ decay is about 1.24 x 
10~^, the direct CP asymmetry Cj^t^ is about -0.05, while the mixing-induced CP asymmetry is 
as large as 5'7r7r=0.18. Our results might be tested in the near future at LHCb. 

The remaining parts of this paper are organized as follows. In the next section, we outline 
the necessary ingredients of the QCD factorization approach for describing the Bs tt+tt^ 
decay and calculate the decay amplitude. In section 3, we give the numerical results of the CP 
averaged branching ratio and discuss CP asymmetries in Bs n^ir" decay. 

2 Bs ^ 7T^7T~ decay in QCD factorization approach 

We will start as usual from the effective Hamiltonian for the AB = 1 transitions given by ^Hj 

--y 10 

He// = ^{K.v;;[C'i(/i)Oi(/.) + c'2(/i)02(^)]-F*,v^,:Ec'.(/^)o.(/.)} + /i.c, (i) 

V-^ i=3 

where Cj are Wilson coefficients at the renormalization scale /i in the Standard Model by 
integrating out heavy gauge bosons and top quark fields. O12 are tree operations arising from 
W-boson exchange and 03_io are penguin operators. The values for Q and the definition of 
operators Oi could be found in [THj. 

With the effective Hamiltonian, the amplitude for Bg tt+tt" in naive factorization is 
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+ ("5 + -Q'aa'j){'n^Tx'\u-i^Ru\Q) 
+ ^^Kfe\/*j5^p^a2(7r+7r-|n7;.Ln|0) + (n ^ d), (2) 



where L,R = (1 =p 75)/2. Due to the conservation of vector current and partial conservation 
of axial-vector current, this amplitude will vanish in the limit mu,md 0. To as order, the 
matrix (tt+tt^Im /' ^(1 — l5)u\0) also vanishes due to the cancellation between the amplitudes 
of Fig.l (a) and (b). So that, nonfactorizable contribution will dominate the decay, which can 
be obtained by calculating the amplitudes of Fig.l (c) and (d). We consider the contribution 
up to the twist-3 distribution amplitude of the light mesons which is superficially suppressed 
by n-ir, however, fi^ is much larger than its naive scaling estimation Aqcd |1] 
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The amplitudes are calculated to be 

G p 



dL 



dx f dy\ 
Jo 



Dskl 



^2 



m 



B 



y<l>^(/+) + 3e$^(/ 



nib 



2— 

rriB 



Dbkl 



(4) 



+ 
+ 



dl+ dx 



2C4 + 



G 



10 



2C6 + ylK^B 



Dskl 



Cs 
2 



2C4 -|- 2Cq 



10 \ /^TT 



07r(a;)0^(?/) 



TUB 



2— ($^(/+) + $^(/+; 



Dskl 



(5) 



where a; = 1 — x, = {Mb — mb)/MB, and ,^ = 1+/Mb. D^^s and /c^ are the virtualities of b 
quark, s quark and gluon propagators respectively. $'s are the leading twist light-cone distribu- 
tion amplitude(DA) of vr and B mesons. (t)T^{x) is the twist-3 DA of tt meson. These distribution 
amplitudes can be found in Refs. [HI [T21 CHI El which describe long-distance QCD dynamics 




of the matrix elements of quarks and mesons, which are factorized out from the perturbative 
short-distance interactions in the hard scatting kernels. For the distribution functions of B 
meson, we use the model proposed in jTT| 
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Now we can write the total decay amplitude 

where z = \ VtbV^*/VubV*J\A^/A^\, 7 = axg[VtbV^*jVubV*g], 6 is the relative strong phase between 
penguin and tree contribution amplitudes, z and 6 can be calculate within QCD factorization 
framework. 



3 Numerical results and Summary 

We list the parameters used in our numerical calculation J3] 

Mb = 5.37 GeV, nib = 4.66 GeV, r^o = 1.461 ps, fs = 236 MeV, 

(9) 

U = 130 MeV, p = 0.20, f] = 0.33. 



We set the scale in as{n) to be Mb J 2 which is about the averaged virtuahty of the time-hke 
gluon. In EqlH El we meet endpoint divergences, which is the known difficulty to deal with 
the annihilation diagram within QCD factorization framework. Instead of the widely used 
treatment Jq y — > Xa = (1 + Qa^^'^) In ^ in the literature [121 El d], we use an effective gluon 
propagator [2] to avoid enhancements in the soft end point region 



1 1 -.n,.n. n rln(^^-4^)l -TT 



F k'^ + MW 



(10) 



Typically rUg = 500 ± 200 MeV, A = Aqcd=250 MeV. The use of this gluon propagator is 
supported by lattice result [1^], and field theoretical studies [3120] which have shown that the 
gluon propagator is not divergent as fast as p-. 

For twist-3 DA 07r(x), its asymptotic form is (f)n{x) = 1 ^B] which used in [Old]. To further 
suppress endpoint contributions, we will use the recent model by Huang and Wu[ 
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where 6*2^^(1 — 2x) and Cy^{l — 2x) are Gegenbauer polynomials and other parameters could 
be found in jT3]. 

Using these inputs, we get the CP averaged branching ratio of the decay 

Br{Bs 7r+7r~) = (1.24 ± 0.28) x 10^^. (12) 

The available upper limit of the decay at 90% confidence level [T3] is 

Br(Bs TT+TT-) < 1.7 X 10"^ (13) 

Obviously, our result is far below this upper limit. However, our result is larger than these 
QCD factorization result Br{Bs tt+tt^) ~ 2 x 10~^ CH! by using the treatment ^ — >• 
Xa = (1 + Qa^^'^) l^^f"' note that our result may consistent with the one of Ref. jH] 

Br{Bs TT+TT^) = (0.024l[]:[J[}^l[j;g?^l[j:J^f) x lO"*' if the huge uncertainties are considered. In a 
recent study in the framework of PQCD factorization [21], the authors found Br{Bs tt+tt^) = 
(4.2 ± 0.6) X 10^^ where the end point divergence is regulated by k"^. 

The absolute ratio between the amplitude of penguin and the tree is 2; = 9.8, and the strong 
phase is 5 = 164°. So, we can see that almost all the contribution comes from penguin. Our 
results for z and 6 agree with the PQCD results |^. 
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Figure 2: The direct CP violation parameter Ctttt and the mixing-induced CP violation parameter 
Stttt of Bg 7r"''7r~ decay as a function of weak phase 7(in degree). 

Now it is time to discuss CP asymmetries of Bs{Bs) ti^h^ . The time-dependent asym- 
metries are given by [22] 

where Am is the mass difference of the two mass eigenstates of meson. C^^^^ and S^y^^ cire pa- 
rameters describing the direct CP violation and the mixing-induced CP violation, respectively. 

Finally, our results for direct and mixing- induced CP violations in the decay are presented 
as functions of weak phase 7 in Fig. 2. a, b respectively. For 7 = 60° ± 14° jl^, the direct CP 
violation parameter Ctttt is about -0.05, the mixing-induced CP violation parameter ^tttt of the 
decay is as large as 0.18. 

In summary, we have calculated the CP averaged branching ratio and CP asymmetries of 
the decay Bs tt+tt^ within the framework of QCD factorization. We have obtained that the 
CP averaged branching ratio of this decay mode is of the order of 10"^. The CP violations 
are estimated to be Ctttt = —0.05, S'^^tt = 0.18. Compared with former studies in the same 
framework, we have included both the two distribution functions and of Bs meson. We 
also have used Cornwall prescription [Zj for the gluon propagator with a dynamical mass to avoid 
enhancements in soft endpoint region. It is noted that recent studies jHl E] have given support 
for Cornwall prescription, which might have many phenomenological applications in B decays. 
Once future measurements at LHCb in agreement with our predictions, it would indicate that 



Cornwall prescription could be used in QCDF to improve it's treatment of endpoint divergences 
in hard-spectator scattering and annihilation topologies to enhance its power for analyzing 
charmless B nonleptonic decays. 
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